Reading-out the state of a flux qubit by Josephson transmission line solitons 
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I. INTRODUCTION 

The last few years have brought substantial break- 
throughs in experiments with Josephson qubits, with fur- 
ther progress depending to a large extent on our ability to 
control the circuits with high precision and, at the same 
time, avoid decoherence. One of the promising ideas com- 
bines Josephson flux qubits with the well developed clas- 
sical rapid single flux quantum (RSFQ) technology 
RSFQ elements should allow for more reliable on-chip 
control and measurements than what can be achieved 
with control pulses sent over long coaxial cables. How- 
ever, they have the drawback to be dissipative and thus 
noisy. The effect of this noise needs to be investigated, 
and ways need to be found to minimize it. 

In this article we consider one of the possible elements 
of an RSFQ circuit, namely a Josephson transmission 
line (JTL). The JTL supports propagating signals in the 
form of Josephson solitons (phase slips), also called flux- 
ons. As suggested by Averin et al. [1] ballistic fluxons 
in the JTL can be used to read out the state of a su- 
perconducting flux qubit in a setup as shown in Fig. [T] 
In the proposed schemes the information about the state 
of a qubit is contained either in the fluxon transmission 
probability (transmission detection mode) or propaga- 
tion time (delay time detection mode). Under certain 
ideal conditions the measurement time is equal to the 
back-action dephasing time, which proves that the JTL 
can operate as an ideal detector. In this paper we investi- 
gate the efficiency that can be achieved in the delay time 
detection mode. For that purpose we evaluate the delay 
time in three different setups: a) when the qubit is kept 
away from the symmetry point all the time; b) when the 
qubit is initially prepared at the symmetry point, but 
the approaching soliton pushes the qubit far from the 
symmetry point; c) when the qubit is near the symmetry 
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FIG. 1: Setup for the read-out of the persistent current qubit 
based on the delay time of a soliton in the Josephson trans- 
mission line (JTL). 



point all the time. We analyze the relation between the 
delay time and the dissipation as well as the probabil- 
ity of errors introduced by the measurement. Finally, we 
compare the delay time with the characteristic time un- 
certainty due to jitter (thermal fluctuations) in the JTL, 
and we determine how many solitons are needed for a 
reliable measurement. 

For a flux qubit operated far from the symmetry point 
the eigenstates are persistent current states. The persis- 
tent current in the qubit loop induces an external mag- 
netic flux in the JTL, which provides a scattering poten- 
tial for the fluxon and is responsible for the time delay 
of the fluxon propagation. The sign of the magnetic flux 
and the value of the delay time depend on the state of 
the qubit, which allows measuring its quantum state. 

For a qubit prepared in one of the energy eigenstates 
at the symmetry point the expectation value of the cur- 
rent in the loop vanishes. However, for strong qubit- JTL 
coupling the fluxon shifts the working point of the qubit 
away from the symmetry point, and then again the qubit 
produces a magnetic flux in the JTL. The delay time is 
approximately the same as if the qubit would have been 
in a persistent current state corresponding to the induced 
shift. 

For weak qubit-JTL coupling a fluxon shifts the work- 
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FIG. 2: Energy diagram of the qubit for the considered mea- 
surement schemes: a) The qubit is kept away from the symme- 
try point, b) The qubit is initially prepared at the symmetry 
point, but the approaching soliton pushes the qubit far from 
the symmetry point, c) The qubit remains near the symmetry 
point all the time. 



ing point of the qubit only slightly around the symmetry 
point. For this case we demonstrate that the qubit in- 
troduces an effective inductance to the JTL with sign 
depending on the qubit eigenstate. The local change of 
the inductance of the JTL also serves as a scattering po- 
tential for the fluxons, leading to a delay in the fluxon 
propagation time. This property can in principal be used 
for a measurement of the qubit at the symmetry point 
even for weak qubit- JTL coupling. 

For a realistic assessment of the feasibility for the pro- 
posed measurement schemes we need to consider the ma- 
jor sources of errors. One of them is thermal noise in 
the JTL, as a result of which fluxons show fluctuations 
in their velocity, which lead to an uncertainty (jitter) in 
the propagation time. In order to distinguish different 
eigenstates of the qubit the induced delay time should 
exceed this time jitter. Another source of errors are non- 
adiabatic transitions of the qubit caused by the moving 
fluxons. This back-action effect of the JTL on the qubit 
needs to be taken into account if the qubit is measured at 
the symmetry point. Still another source of errors are in- 
trinsic relaxation processes of the qubit due to noise not 
related to the JTL. The relaxation destroys the state of 
the qubit to be measured and should be much slower than 
the measurement time. Finally, in order to be detectable 
the delay times should exceed the time resolution of the 
RSFQ detector used for delay time measurement. 

Our calculations show that for suitable JTL parame- 
ters the induced delay times well exceed the time reso- 
lution of the RSFQ detector. For strong qubit- JTL cou- 
pling the qubit measurement by a single fluxon at the 
symmetry point, as well as far from it, can be performed 
with an accuracy of 70 — 90% for experimentally accessi- 
ble parameters. For high velocities of the fluxon the mea- 
surement efficiency is mostly limited by the jitter and, if 
the qubit is at the symmetry point, by Landau-Zener 
transitions. For low velocities the intrinsic relaxation of 
the qubit may be important. To improve the accuracy 



one has to use many fluxons for the measurement. If one 
optimizes the fluxon velocity and number the qubit can 
be measured with accuracy reaching 99% far from the 
symmetry point and above 90% at the symmetry point 
for strong qubit- JTL coupling. The measurement time 
remains approximately 4 ns for both cases. For weak 
qubit- JTL couping the qubit prepared at the symmetry 
point stays there during the whole measurement. In this 
case the measurement of the qubit is not feasible for stan- 
dard JTL parameters. 



II. THE MODEL 

We focus on the readout of a persistent current 
qubit I, 3 via its coupling to solitons in a Josephson 
transmission line (JTL) in a setup as shown in Fig. [1] 
The qubit consists of a superconducting loop with three 
Josephson junctions and is inductively coupled to the 
JTL. The Hamiltonian of the total system is 



H = H. 
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where Hqb describes the qubit, Hjtl the JTL, and Hj 
the qubit-JTL coupling. 

In the basis of the two lowest persistent current states, 
|0), |1), corresponding to currents ±Ip circulating in the 
qubit loop in opposite directions, the qubit Hamiltonian 
can be expressed in terms of Pauli matrices Q , 
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The tunneling amplitude A, leading to transitions be- 
tween both states, is fixed by the experimental setup. 
The energy bias Cq between the two persistent current 
states depends on the deviation of the external magnetic 
flux <f> in the loop from the symmetry point. 



eo = 2/p($-$o/2), 



(3) 



where $o is the magnetic flux quantum. 

The Hamiltonian of the discrete dissipationless JTL, 
shown in Fig. [31 in vanishing external magnetic field can 
be expressed in terms of the charges g„ and the phase 
differences 0„ of the n**^ Josephson junctions 
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Here, N is the total number of junctions in the JTL, 
which here are assumed to be equal, C is the capacitance, 
and Ej — hlc/2e the Josephson energy of the junction 
with Ic being the junction critical current. The induc- 
tance of each cell of the JTL is denoted by L, and /g is 
the bias current supplied externally to each junction. 
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The dimensions of the qubit are much smaller than the 
length a of a unit cell of the JTL. Therefore, one can as- 
sume that the qubit is inductively connected only to one 
cell of the JTL denoted by label m. The corresponding 
mutual inductance is M — kyjLqhL, where Lqh is the 
inductance of the qubit loop and k is the coefhcient of 
coupling. Thus the Hamiltonian of the qubit-JTL inter- 
action in the linear response regime is given by 



Hi = a.MIpiL 
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where is the electrical current in the inductance of the 
m-th JTL cell and g = 2JpM$o/(27ri) is the effective 
qubit-JTL coupling strength. 

The classical dynamics of the uniform, discrete, dissi- 
pative JTL is governed by Kirchhoff circuit equations [5] 
which lead to a system of N discrete sine-Gordon equa- 
tions 
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Here, 



(2e/c/fi.C)"'^/^ is the plasma frequency of the 



junction, lOc = 2eIc/hGN the junction characteristic fre- 
quency with Gn being the normal conductance of the 
junction determining the dissipation in the system, and 
<i>^ is the external magnetic flvix in the n— th cell induced 
by the qubit. 

For low inductances. Lie ^Oi and weak external 
magnetic fields, $^ <C $o, the phases <j)n of the neighbor- 
ing junctions are close to each other, and we can rewrite 
([^ in dimensionless differential form. 



(7) 



Here the time t and continuous space variable x 
are measured in the units of the inverse plasma fre- 
quency cj"^ and the Josephson penetration depth Aj = 

1/2 

a ($o/(27rL Jc)) , respectively, and the phase difference 
<j){x) is a function of x. Furthermore, a = LOpjiOc is the 
dissipation strength due to a non- vanishing normal con- 
ductance Gat, and f'^''{(p,x) is the perturbation induced 
by the qubit. 

In the following sections we will show that far from the 
symmetry point, Eq 3> A, the perturbation is created by 
the magnetic flux due to the qubit and can be expressed 
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where the different signs, ±, correspond to the two persis- 
tent current states |0) and |1) of the qubit, respectively. 

For a qubit remaining at the symmetry point, eo = 0, 
g ^ A, we obtain an inductive type interaction, and the 
perturbation is given by 



LXjja dx 
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FIG. 3: Two elementary cells of the discrete JTL. 



Herethe different signs, ±, correspond to the energy 
eigenstates of the qubit at the symmetry point, |±) = 
(1/V2)(|0)±|1)). 

In the limit a = j e = f^^ — the exact soliton solution 
of Eq. ([7]) has the form 



4>o{x, t) = 4 tan 
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where the positive (negative) sign corresponds to a fluxon 
(antifluxon), is the initial position of the soliton, and 
u its velocity, which can take any value between —1 and 
1. To analyze the fluxons dynamics in the JTL follow- 
ing from Eq. ([7]) in the general case we make use of the 
collective coordinate perturbation theory developed by 
McLaughlin and Scott @. It is based on the assump- 
tion that je, ct.(j), x) <C 1, which allows writing the 
fluxon solution in the form 



4>q{x, t; u, X) — 4 tan 
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For simplicity we consider here the case where only one 
fluxon is present in the system. 

Without coupling to the qubit, = 0, the stationary 
velocity of a fluxon can be derived from the power balance 
equation [gi], and is given by 
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When the qubit is coupled to the JTL, for arbitrary form 
of the perturbation /'^''(0, x), the variation parameters 
X{t) and u{t) obey the differential equations [6| 
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where = (a;— X(i))/(1— ■u^)^/^. In the following section 
Eqs. and will be solved for both forms of the 
qubit perturbation /'^ given in Eqs. ([5]) or ([5]). 
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TABLE I: Parameters of the JTL and the persistent current qubit. 





Parameters of JTL [1, 10] 


j Qubit-JTL coupling, k 


Parameters of persistent current qubit [12] 


/c(mA) 


MA/cm^) Xj{a) L(pH) cjp/27r(GHz) 


1 »>A g<A 


A/27r(GHz) e™„,/27r(GHz) L,6(pH) /j,(nA) 


0.6 


30 2 137 24 


i 0.96 0.05 


5.5 55 5 300 


2 


30 2 41 24 


1 0.53 0.025 





III. DELAY TIMES. 



dimensionless coupling 



A. Qubit far from the Symmetry Point 

We first consider the situation where the qubit is pre- 
pared far from the symmetry point, S> A. In this 
case the eigenstates of the qubit are the persistent cur- 
rent states, |0) and |1). For each of them a magnetic flux 
zLMIp penetrates the cell of the JTL which is inductively 
coupled to the qubit. If the size of the qubit is much less 
than the Josephson penetration depth Aj the external 
magnetic flux in the JTL can be written as 



< 



$^(a;) = ±MIpe{a 



(15) 



where 9{x) is the step function, and we assumed the qubit 
to be located at a; = 0. The corresponding perturbation 
term in the equation of motion ([7]) can be written as 



(16) 



where (j>l = 27rM/p/$o is the dimensionless qubit-JTL 
coupling. A thorough experimental and theoretical study 
of the fluxon dynamics described by the sine-Gordon 
equation with delta-function terms was performed in 
Refs. ff', '8*1 for long annular Josephson junctions. Propa- 
gation of the charge soliton in ID array of serially coupled 
Josephson junction was studied in Ref. @. 

The dimensionless coupling (pf depends linearly on the 
mutual inductance M, which for typical experimental 
conditions is much lower than 1. Increasing the cou- 
pling coefficient k between the qubit loop and the JTL 
cell enhances the influence of the qubit on the fluxon dy- 
namics in the JTL, but at the same time increases the 
back-action of the fluxon on the qubit. In particular, the 
fluxon induces a magnetic flux in the qubit loop, which 
shifts the working point of the qubit. The persistent cur- 
rent qubit is only well defined when the external flux 
in the qubit loop is close to the symmetry point $o/2. 
For larger values of the external magnetic flux the form of 
qubit-JTL perturbation (|16p is not valid and a more elab- 
orated model should be used. Moreover, large deviations 
of the working point can lead to non-adiabatic transitions 
to higher levels, which also makes the two-level model in- 
valid. In the following we require Mlpl^lmax < emaa;/2, 
where emax is the maximum value of the energy bias be- 
tween two persistent current states for which our system 
is still regarded as the persistent current qubit. This 
condition translates to the following limitation for the 
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and for the coupling coefficient k — M/ y/LLgt, 
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The right hand side of Eq. ^T7\ is the ratio between the 
qubit energy splitting and the magnetic energy of a fluxon 
<I'oa/(7r^LAj). It demonstrates that the efficiency of the 
measurement is low if the magnetic energy of a soliton is 
much higher than the qubit energy, because one needs to 
decrease the coupling coefficient k according to p^ . 

The magnetic energy of a soliton in the JTL can be 
reduced by decreasing the critical current Ic of the junc- 
tions. We evaluate the delay times for two specific val- 
ues of the critical current, Ic = 0.6 ^A and Ic = 2 /iA. 
Given these values the cell inductance L has to be cho- 
sen appropriate to yield suitable values of the Josephson 
penetration depth Aj. This length should be large to 
minimize effects of the discreteness of the JTL, Aj > 2a. 
On the other hand, it is useful to have Aj as small as 
possible to decrease the magnetic energy of a fluxon and 
to keep the number of Josephson junctions in the JTL 
within practical limits. Based on these considerations [l| 
we choose Aj = 2a. The parameters of the JTL and the 
qubit which were used in the calculations are listed in 
Table D 

From (|13m6p we obtain 



±i(l - u2)i/20«sech^eosheo, 
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dX 



u ± -u(j)l (coshGo - eoshGo) sech^Go. (20) 



where Gq = X/{l-u^Y'^. 

In the ballistic regime, je = and a = 0, we can inte- 
grate (froi I^O)) analytically and obtain the fluxon velocity 
u and center coordinate X as functions of the parameter 
00 for the qubit prepared in the eigenstates |0): 
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The solution for the state |1) is obtained by replacing (f>1 
by -</>?. 

To proceed we first consider the situation, where the 
soliton is able to pass the potential and, therefore, we 
can always introduce the delay time caused by the per- 
turbation due to the qubit as 

td = |^|o) ~ = 2|t|o) - ^ol 



2 r°° 




UO 7-00 


u{X) 



dX. (23) 



Here t\Q) and t\i) are the propagation times correspond- 
ing to the two persistent current states of the qubit, and 
to is the value without qubit interaction. The factor 2 in 
Eq. ([^ reflects the fact that the perturbation /g^^ has 
the same magnitude but different signs for the two per- 
sistent current states. When the velocity is only slightly 
perturbed by the qubit, ei = (1 — Uq)4>\/ {Auq) <^ 1, we 
find 

u(eo) =Mo(l-£isecheo)-t-o(£?). (24) 
The delay time ((23|) can then be evaluated as 
2s / \ 

^ ^ loo ^^"'^^ I (T^^ j " 4' ^''^ 

where the last part of the equation describes the "non- 
relativistic case" uq <^\. 

From (HH) we see that for (f)\/A > (1 - ul)-^/'^ - 1 the 
soliton does not have enough kinetic energy to pass the 
potential barrier induced by the the qubit in the state |0). 
After approaching the qubit the soliton will be reflected 
by the barrier. Due to the combination of dissipation and 
driving the soliton will oscillate and eventually stop at 
the "pinning point" where u = Q and X = Q. The pinning 
of a fluxon and can also be used for the transmission 
detection mode of measurement, but will not be analyzed 
here further. 

The results of a numerical evaluation of the delay times 
for the JTL with dissipation are presented in Fig. 01 As 
expected the analytical formula (pS)) gives a good esti- 
mate even for intermediate values of the soliton velocity 
Wp ^ mq <C 1, where Up is the minimal possible veloc- 
ity for which a fluxon is still able to pass the potential 
barrier created by the qubit. 

Another observation is, that for fixed value of the soli- 
ton velocity uq the delay time decreases with increasing 
strength of the dissipation a. This property is displayed 
by Fig. [51 which shows the velocity as a function of the 
soliton coordinate for different values of a. One can see 
that with increasing dissipation the fluxon velocity devi- 
ates from the dissipationless solution (|13ll4p . The initial 
decline of the fluxon velocity becomes less effective and is 
also partially compensated by the following acceleration 
of the fluxon. The compensation is the more complete 
the higher is the dissipation, which leads to zero delay 
times in the limit of strong dissipation. This behavior is 
typical for the particle driven in viscous media and can 



be understood from power balance considerations. Let us 
consider the case when the qubit induces a positive po- 
tential barrier for a fluxon. When a fluxon is moving with 
stationary velocity far from the qubit its energy gain and 
losses are equal. Approaching the potential barrier the 
velocity of a fluxon decreases, which leads to a reduction 
in energy losses due to dissipation. However, the gain of 
energy due to driving stays on the same level. As a re- 
sult a fluxon receives an excess amount of energy, which 
later leads to an increase of the fluxon velocity above the 
initial value. Thus, we can use the dissipationless solu- 
tion as an approximation for evaluating the delay times 
of a fluxon only if a < 0^. For a > the delay times 
degrade rapidly with increasing a. We expect that the 
delay times can be reliably detected only if a < 0^ . This 
establishes another limitation to the parameters of the 
JTL for the case of very weak qubit- JTL coupling. 



Next we consider the situation when the qubit is ini- 
tially prepared at the symmetry point but shifted far 
from it by the moving fluxon. In what follows we will 
show that this regime is approximately equivalent to that 
when the qubit is kept far from the symmetry point all 
the time, and we can use the previous results for the 
delay time. Indeed for typical parameters of the JTL 
and chosen inductive coupling we have 5 3> A, which 
proves that even if the qubit is initially prepared at 
the symmetry point, eo = 0, it is shifted far from it 
by the flux induced in the qubit loop by the moving 
fluxon. At the symmetry point the persistent current 
in the qubit loop is zero, and the qubit induces no mag- 
netic flux in the JTL cell, lL{t) — Q. As soon as the 
approaching fluxon shifts the qubit from the symmetry 
point, 2MIpI]^{t) ~ A, the persistent current is no longer 
zero. The maximum value of the persistent current, /p, is 
reached once the qubit is pushed far away from the sym- 
metry point, 2IpMlL{t) ^ A. When the fluxon moves 
away the qubit returns to the symmetry point with zero 
persistent current in its loop. These considerations, of 
course, are valid only if the qubit was prepared in an en- 
ergy eigenstate and if the shift by the fluxon occurs adi- 
abatically. The qubit affects the fluxon dynamics only 
during the period of time ~ Xj/uq when the fluxon is 
moving close to the qubit. Most of this time the qubit 
is already shifted far from the symmetry point and in- 
duces the perturbation to the fluxon motion according 
to Eq. ([5]). The quick switching of the persistent current 
from zero to Ip and back in the initial and final stages of 
the qubit- JTL magnetic flux exchange yield only a small 
contribution to the overall delay time experienced by the 
fluxon. Therefore, we can use the calculated delay times 
shown in Fig. 0] even if the qubit is initially prepared at 
the symmetry point for 5 3> A. However, we need to 
evaluate the probability of non-adiabatic transitions of 
the qubit to another eigenstate, which creates an addi- 
tional source of error for the measurement. These results 
will be presented in Section [TVl 
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FIG. 4: Delay time of a soliton induced by the qubit as a function of initial soliton velocity no when the qubit is far from the 
symmetry point, eo ^ A. The velocity and time are measured in units of XjcOp and ujp^ , respectively. Four curves are obtained 
numerically from Eqs. (|13l I14|l for different values of the dissipation strength a (shown in the legend box). The red solid line 
shows the analytical solution (|25|l . The parameters of the JTL and the persistent current qubit are listed in Table U] The inset 
shows the dependence of the delay time on the dissipation strength a for the initial soliton velocity uq = 0.2. 
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fluxon coordinate x, dimensionless units 
FIG. 5: Fluxon velocity as function of its coordinate for differ- 
ent values of the dissipation strength a (shown in the legend 
box) 



B. Qubit at the Symmetry Point 

In this subsection we consider the situation where the 
qubit is prepared initially at the symmetry point eo — 
and remains nearby during the measurement. This re- 
quires that the coupling between qubit and JTL is weak 
and g <C A to ensure that a soliton does not shift the 
qubit from the symmetry point significantly. As a con- 
sequence the expectation value of the flux in the energy 
eigenstates of the qubit is close to zero. If the coupling 
term ([S|) varies slowly we can treat the sum Hqi) -\- H I in 



an adiabatic approximation and diagonalize it to obtain 
AyA + M^I^Ilit) (26) 

A 



qb 



H, = 



Pz 



Leffllit) 



Tjadiab 



diab 



jjadia 



where ie// — 2/pAf^/A. The approximation is valid if 
P ^ 1, where P is the probability of Landau- Zener tran- 
sition between energy eigenstates due to time dependence 
of lL{t), which is calculated in the next section. 

The interaction term Hf'^'^"'^ indicates that the qubit- 
JTL interaction in the adiabatic approximation is equiva- 
lent to the change of the inductance of the JTL cell which 
is coupled to the qubit by the additional value of L^ff 
whose sign depends on the energy eigenstate of the qubit. 
This property can be used for the read out of the qubit 
when it is always kept at the symmetry point (f> = $o/2. 

The qubit-JTL interaction Hf'^'"'' leads to the follow- 
ing perturbation term in the sine-Gordon equation ^ 

,^ [<5(x)0,] , (27) 

where (f)\ — L'^^ ^ a/ {L\j) is the corresponding dimen- 
sionless qubit-JTL coupling and ± correspond to the 
eigenstates |±) = (l/^/2)(|0) ± |1)). 

The Hamiltonian (pS)) is valid only for g <C A, which 
establishes the condition for the maximum possible di- 
mensionless coupling 

A/2 



/f„>,a;)=±0l£-[5(x)0,] 



and coupling coefficient 



k < 



AAj 



^qb 



1/2 



(28) 



(29) 
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FIG. 6: Delay time of a soliton induced by the qubit as function of the initial soliton velocity mq when the qubit remains at the 
symmetry point eo = 0. The velocity and time are measured in units of XjtOp and lo^^ , respectively. Four curves are numerical 
solutions of Eqs. (|13l I14p for different values of the dissipation strength a (shown in the legend box). The red solid curve 
shows the analytical solution (|34|l . The parameters of the JTL and persistent current qubit are given in Table U The coupling 
coefficient between the qubit loop and the cell of the JTL is A; = 0.05 which leads to the dimensionless coupling (f)\ = Z x 10^^ 
in the sine-Gordon equation. The value of the coupling coefficient is reduced in order to keep the qubit all the time at the 
symmetry working point. 



In order to keep the qubit at the symmetry point for 
the chosen parameters of the qubit and the JTL (shown 
in Table |T| the coupling coefficient should satisfy: k <C 
0.096 for Ic = 0.6 fik and k < 0.053 for I,, ^ 2.0 a*A. For 
the calculation we take k — 0.05 and k = 0.025 which 
leads to 02 = 3 X 10""* and 4)\ = % x 10"^, respectively. 

From I|13ll4p and (H?]) we obtain 



du 
It 



dt 



± sech^ 6o tanhOo , 



(30) 



u± 



bl u sech^9o 



(1_ ^2)1/2 (l-Q"tanh9o). (31) 



For a = j'e = in the non-relativistic limit w <C 1 one 
can evaluate for the qubit in the eigenstate 1+) 



and 



=ul~ ^sech^X, 



u ~ uo (l — e2sech^X) + o{e\), 



(32) 



(33) 



where £2 = (tA/i'^ul). The pinning of the fluxon occurs 
if Mq < 02/2- We can calculate the delay time according 
to (IH as 



= ^ / sech^XdX = ^. 



Mo J- 



(34) 



The results of the numerical evaluation of the delay 
times in this regime are shown in Fig. [51 For the same 



values of the fluxon speed the delay times are smaller 
as compared to Fig. [J] However, since the pinning of a 
soliton occurs at slower velocities we can, in principle, 
achieve larger values of the delay times. 



IV. BACK-ACTION OF A FLUXON ON THE 
QUBIT 

The back-action of a fluxon on the qubit produces sev- 
eral effects including dephasing of the qubit in the mea- 
surement basis [l^l and a shift of the working point of the 
qubit. Dephasing, i.e., the decay of the off-diagonal ele- 
ments of the density matrix of the qubit in the measure- 
ment basis, does not affect the measurement outcomes 
studied here [T3| and is not considered further. 

In contrast, the shift of the working point of the qubit 
can potentially "destroy" the qubit and, furthermore, in- 
duce non-adiabatic transitions between the qubit states. 
Both processes create errors in the measurement and 
should be taken into account. By restricting the qubit- 
JTL coupling according to (fTSjl or (P^)) we insure that 
the system remains a persistent current qubit during the 
whole time of the measurement. In this section we derive 
the probability of non-adiabatic transition between eigen- 
states of the qubit due to a passing fluxon. It should be 
noted that non-adiabatic transitions are a potential prob- 
lem only if the qubit is initially prepared at the symmetry 
point. Far from the symmetry point the Hamiltonian of 
the qubit ^ and the interaction Hamiltonian ^ ap- 
proximately commute and non-adiabatic transitions are 
negligible. 
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In order to study nonadiabatic transitions we assume 
that the qubit is initiaUy prepared in the excited state 
|— ) at the symmetry point eo — where its unper- 
turbed dynamics is described by the Hamiltonian Hqb = 
— {A/2)ax- For a fluxon passing the qubit the time- 
dependent qubit- JTL couphng can be written according 
to IS]) and Uni) as 



Hi = MIplL{t)c7, 



(35) 



■nL{\j/a)^l - ul 



sech 



Here we neglected the action of the qubit on the fluxon. 
The transition probabihty from the excited state |— ) to 
the ground state |-f) induced by one passing fluxon is 
then given by 



p = Mip{2hy 



sin A 



A 



(36) 



where A = MIp{2h)-'^ J^lL{t)dt. After inserting ^ 
we obtain 



sm 



L(Xj / a)uohujp 



sech 



(37) 

We wifl use this result to evaluate the measurement effi.- 
ciency in Section IVlII 

V. SOLITON JITTER 

In order to be detected the delay time should exceed 
the time fluctuations induced by thermal noise in the 
JTL. We can account for thermal noise by adding a 
stochastic term ^{x,t) to the sine-Gordon equation ([7]), 



(38) 



The noise is assumed to be white with autocorrelation 
function 

{^{x, t)^{x',t')) = 16a{kBT/Eo)d{x - x')5{t - t'). (39) 

Here i?o = 8{hJcXj/2e) is the rest energy of a soli- 
ton, and the coefflcient IQajh sT l Ep) is fixed by the 
fluctuation-dissipation theorem [iq . The noise spectral 
density for the fluxon velocity is given by [l^ 



2akBT{l-u^) 



2^5/2 



Eo 



(40) 



where Au(t) = u{t) — uq. Starting from the initial condi- 
tion that at i = the velocity of the fluxon is controlled 
and equal to uq we find the velocity autocorrelation func- 
tion to be given by 



{Au{h)Au{t2)) 



-(1 



^g-Q|tl-t2| _ g-a|ti+t2| 



.(41) 



The uncertainty in the velocity of the fluxon leads to 
an uncertainty of its coordinate AX{t) = J* Au{t')dt' 
and propagation time At — AX/uq. The time jitter 
of the fluxon, defined as the standard deviation of the 

1/2 

propagation time St = ((Ai^)) , is then given by 



St 



X \ta + e 



Uo 



1/2 



(42) 



Limiting cases of Eq. are St oc t^^^ in the diffusive 
regime 3> 1, and St oc t^^'^ in the ballistic regime 
at ^ 1. The time jitter (I42p leads to further errors in 
the measurement procedure, which will be analyzed in 
section VII. 



VI. RSFQ DELAY DETECTOR 

Before turning to a quantitative analysis of the mea- 
surement errors we describe in this section an RSFQ de- 
lay detector that is need to measures the time between 
two SFQ pulses propagating through a JTL. We first eval- 
uate the time resolution of a single RSFQ decision gate 
and the time resolution of the improved detector based 
on a time vernier. Based on realistic parameters we es- 
timate the hardware complexity of the detector and dis- 
cuss possible designs of the circuit presently developed 
experimentally [19|. 

The operation of any RSFQ gate is based on the time 
resolution between SFQ pulses. The gates produce bi- 
nary output depending on the relative time between ei- 
ther clock and data pulses or between two input data 
pulses. The simplest RSFQ gate that can be used as 
a decision circuit is the asynchronous OR or confluence 
buffer shown in Fig. [T] The flgure also illustrates the 
operation of the confluence buffer as a time detector. A 
single SFQ pulse on either input produces a SFQ pulse on 
output. Therefore, for delayed input pulses two output 
SFQ pulses are produced. In case of simultaneous arrival 
the input pulses compensate each other resulting in a 
single SFQ pulse on output. Such behavior of the RSFQ 
confluence buffer is well known and has been experimen- 
tally verified many times, see for example Ref. |20| |. 

In general, the time resolution of any RSFQ gate is 
determined by the resolution of the balanced compara- 
tor formed by two junctions connected in series. In the 
confluence buffer balanced comparators are formed by 
junctions Jl, J3 or J2, J3. The probability of junction 
switching in the comparator obeys a normal distribution 
with CTc ~ O.lStsFQ, where tsFQ is the width of SFQ 
pulse [21], 123]. In time units normalized to 1/ujp, the 
width of the SFQ pulse is approximately 



tsFQ ~ 6- 



1 



(43) 
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FIG. 7: (a) Schematic of RSFQ confluence buffer, and (b) the probabihty of double switching of output junction J3. 



where = ^^c/^p is the McCumber-Stewart parame- 
ter. The corresponding normahzed one-sigma jitter of a 
single decision gate with /3c = 2 is 



Oc = G.msFQ ~ 0.78- 



0.55. 



(44) 



For a measurement error of 10~^ the time resolution of 
a single RSFQ gate, t^, is given by 4ctc and equal to 

— 2.2. This time resolution does not dependent on 
the temperature since the measurement involves symmet- 
ric stochastic processes of the switching of two balanced 
junctions 23]. 

As it is shown in Section [1111 the delay time between 
SFQ pulses that needs to be detected is in the range of 

= 0.1 — 5 depending on the actual dissipation in long 
JTLs and the accuracy in fixing the initial fluxon veloc- 
ity. In order to improve the time resolution of the RSFQ 
delay detector time a vernier with N decision gates can 
be used 0, , the block diagram of which is shown in 
Fig. [H After the JTLs, both pulses propagate through 
the chain of splitters and arrive on a chain of N decision 
gates with inputs delayed in time by td- The operation 
principle of the time vernier is the same as that of a multi- 
bit analog-to-digital converter where each bit compares 
the signal with a slightly shifted threshold. In accordance 
with the standard theory for analog-to-digital converters 
the resolution improves as the square root of the number 
of bits [1^. 

Time resolution of the vernier is — t!^/VN and 
relative time difference between stages of the vernier is 
ts — t-^/N. The error of time measurements depends on 
the ratio between time difference, ts, and jitter accumu- 
lated in the pulse pathos, <7y, 



ierfc( — ). 



The accumulated jitter is 



(45) 



(46) 



where Ny is number of Josephson junctions in each pulse 
path and aj is a jitter per Josephson junction. The factor 



Reference pulse 




FIG. 8: Block diagram of a RSFQ time vernier. 



two comes from the fact that there are two chains of 
splitters involved. The number of Josephson junctions in 
the splitter chain grows like Ny = NlogN. The single 
junction jitter at 4.2 K is about aj = O.OlStsFQ [HI and 
it scales as square root of temperature and as square root 
of shunts resistors of the junction. 

Typical operating bath temperatures for the qubit ex- 
periments are about 30 mK. For process with operat- 
ing frequency below 5 GHz and assuming use of cool- 
ing fins for thermalization of the hot electrons in the 
normal metal shunts, the effective noise temperature of 
the RSFQ gates can be reduced down to 80 mK for a 
30A/cm^ process and to 30 mK for a lOA/cm^ process 
pTt - The temperature range between 30-80 mK corre- 
sponds to one-sigma jitter of crj — 0.0012 — 0.0004 for a 
single Josephson junction with /3c = 2. 

For measurement errors of 10~^ the time resolution of 
the RSFQ time vernier is given by Aay that results in the 
following relation for the optimum number of stages. 



^ = ^ N log Na J. 



(47) 



Equation (|T7)) gives = 40 and corresponding = 0.33 
for 30mK (10 A/cm^ process) and A^ = 20 and corre- 
sponding = 0.5 for 80 mK (30 A/cm^ process). 
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VII. MEASUREMENT ERRORS 



25]), (IMI) and as 



In an experimental implementation the length Ljtl 
of the JTL is finite. For Ljtl ^ Aj the delay time 
of a fluxon induced by the qubit is that of the infi- 
nite JTL shown in Figs. [4] and [6l The time jitter of a 
fluxon at the end of the JTL is given by Eq. (j42|) with 
t — Ljtl/uq. Due to the jitter the propagation time of 
a soliton will be scattered around the average values i|o) 
and corresponding to the different energy eigenstates 
of the qubit |0) and |1) if the qubit is far from the sym- 
metry point. (A similar analysis holds for distinguishing 
the states |-t-) and |— ) if the qubit is at the symmetry 
point.) We can introduce the corresponding distribution 
functions P\o)it) and P\i){t) for the propagation time of 
a soliton, which have maxima at t|o) and t^i^, separated 
by the delay time td and whose standard deviations are 
equal to the jitter St. The overlap of the two distribution 
functions P\o){t) and P\i){t) characterizes the potential 
to distinguish the two eigenstates of the qubit. For a 
normal distribution the area of the overlap and the as- 
sociated error of the measurement are eTic[td/ {2V26t)] 
and (l/2)erfc[td/ (2\/2(5i)], respectively For experimen- 
tally relevant case of large separation of the distribution 
functions, t^ ^ 6t, the error can be approximated by 
[V2St/td) exp {-tySSt'^). The ratio of the delay time to 
the jitter, which determines the error of measurement due 
to thermal fluctuations, can be estimated for a JTL with 
low dissipation, a <C (1)1 2 ^ according to the expressions 



td 

(x 

5t 2ul'^ \2kBTLjTL 



1/2 



(48) 
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FIG. 9: Sketch of the distribution functions which describe 
the propagation time of a soliton. The two functions P|o) (f) 
and P|i) (t) are peaked around the average values of the propa- 
gation times i|o>,|i) corresponding to the different eigenstates 
of the qubit |0) and The standard deviation and peak sep- 
aration are given by ^JnSt and ntd, respectively. The overlap 
of the distribution functions (grey color) indicates the regime 
where the measurement yields an error. 



This results suggests that in order to decrease the effect 
of thermal noise it is favorable to use slower fluxons. On 
the other hand, low velocity fluxons can be difficult to 
produce and control technologically. To avoid pinning by 
the qubit we also cannot choose the speed of the fluxons 
too small. 

For realistic experimental situations some of the mag- 
netic flux from the control circuit of the qubit '^ext pene- 
trates the JTL and create an additional potential barrier 
for fluxons. If the external flux '^ext is much smaller 
than the flux quantum <I>o we can again use the collec- 
tive perturbation theory to analyze its effect on the delay 
time. According to the definition ([23)1 there is no effect 
of $ea;t ou the difference in the delay times since the po- 
tential barrier does not depend on the qubit quantum 
state. However, for ^ext > MIp the pinning velocity of a 
fluxon will be larger. This should be taken into account 
and can affect the measurement because low fluxon ve- 
locities could be no longer accessible. 

For further improvement of the signal-to-noise ratio 
one can use n fluxons, sending them to the JTL one af- 
ter another. After n fluxon passings the sum of the de- 
lay times and separation of the peaks of the distribution 
functions will scale linearly as ntd, but the jitter will 
increase only as y/nSt, and the resulting measurement 
error decreases according to {l/2)eik[y/ntd/ {2^/26t)]. 

The qubit also experiences other sources of noise dur- 
ing the measurement which are not related to the JTL. 
Their effect can be described phenomenologically by two 
characteristic time scales: the relaxation time Ti and the 
decoherence time T2. The pure decoherence, i.e., the de- 
cay of the off-diagonal elements of the density matrix 
of the qubit, has no effect on the measurement outcome 
and is not relevant in our case. The relaxation changes 
the probability distribution described by the diagonal el- 
ements of the density matrix and, thus, affects the mea- 
surement. For a qubit initially prepared in the excited 
state, after n passings fluxons the probability to find the 
qubit in the ground state is 1 — exp{—nL / {uqTi)) . We 
take this probability as the approximate (upper bound) 
measure of the possible error of measurement due to re- 
laxation of the qubit. 

Finally, if measure the qubit at the symmetry point 
there is an additional error associated with Landau-Zener 
transitions between the qubit eigenstates. Following con- 
siderations similar to the previous paragraph we define 
the measure of the error due to Landau-Zener transitions 
as 1 — exp(— nP) where probability P is given by ([57)1 . 

The efficiency of the qubit measurement by a single 
fluxon is given in Table HIl The qubit can be measured at 
the symmetry point and far from it with accuracy of ap- 
proximately 90% for Ic = 0.6 fiA and 75% for Ic = 2 fiA. 
To improve the accuracy one needs to use many flux- 
ons. Figs. [TUl and [TT] show the error of the measurement 
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FIG. 10: Different contributions to the measurement error for Lr — 25Aj, T — 20 mK, a — 0.001. The upper plots correspond 
to 7c = 0.6 and k — 0.96 while for the lower plots Ic — 2 /lA and k — 0.53. The qubit is initially prepared far from the 
symmetry point. See text for detailed explanations. 



for the qubit initially prepared far from symmetry point 
and at the symmetry point, respectively. The total er- 
rors of the measurements are shown on the right contour 
plots as a function of the number of fluxons passings and 
the initial velocity of the fluxons. They are estimated 



TABLE II: Error of measurement 


by a single 


fluxon for 


LjTL = 25A,7, T = 20 mK, uo = 0.2, 


a = 0.001. 




Variant 1 


Variant 2 


Variant 3 


/c(mA) 0.6 


2 


0.6 


k 0.93 


0.53 


0.05 


td(t^p') 4.5 


1.2 


0.04 


5t{Lj-^) 1.5 


0.84 


1.5 


Errors of measurement (in %) 


jitter 7 


24 


50 


relaxation 0.08 


0.08 


0.08 


LZ transitions, eo = 2 


2 


6.8 


total error, eo 2> A 7 


24 




total error, eo = 9 


26 


55 



as 1 - exp [-erfc(V^t<i/2\/2(5t) - nL/{uoTi)] for Fig.[TO] 
and 1 - exp [-eTk{^/ntd/2V2dt) - nL/{uoTi) - nP] for 
Fig. [TT] The different contributions to the total error are 
also plotted separately as a function of the fluxon velocity 
(the left plots) and number of fluxon passings (the insets 
of the left plots). The non-adiabatic transitions are sup- 
pressed if the qubit is prepared far from the symmetry 
point, hence there is no curve related to Landau-Zener 
transitions in Fig. [TOl One can see that for high fluxon 
velocities the measurement quality is limited by the jitter 
and the Landau-Zener transition probability (if the qubit 
is at the symmetry point). For low velocities the mea- 
surement time increases and the relaxation of the qubit 
becomes also important. Under optimum conditions the 
qubit can be measured with probability of error below 1% 
far from symmetry point and below 7% at the symmetry 
point with the measurement time of approximately 4 ns. 

Fig. [T^ shows the measurement error for Ljtl — 25Aj, 
T = 20 mK, a = 0.001 and k = 0.05. The qubit stays at 
the symmetry point during complete measurement. The 
coupling coefficient is chosen small to prevent a shift of 
the qubit from the symmetry point by the passing soli- 
ton. Thus, the influence of the qubit on a soliton is also 
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FIG. 11: Different contributions to the measurement error for Ljtl = 25Aj, T — 20 mK, a — 0.001. The upper plots 
correspond to 7c = 0.6 fiA and k — 0.96 while for the lower plots Ic — 2 fiA and k — 0.53. The qubit is initially prepared at 
the symmetry point but is shifted far from it by a moving soliton. See text for detailed explanations. 



weak, which affects the sensitivity of the measurement. 
From Fig. [T^] it becomes clear that in order to achieve 
better measurements we need low soliton speed and a 
large number of fluxons passings. As a result the time 
of measurement which is required to extract the infor- 
mation about the qubit quantum state is comparable to 
the relaxation time of the qubit. This makes it impossi- 
ble to efficiently measure a qubit in this regime for the 
chosen parameters of the JTL. A reduction of the critical 
current of the Josephson junctions may allow decreasing 
the magnetic energy of the solitons and consequently in- 
crease the coupling coefhcient. In this case the accuracy 
of the measurement will be improved while the qubit will 
remain all the time at the symmetry point. 



n=200 



- thermal fluctuations 
Landau-Zener transitions 
qubit relaxation 




I =0.6 



0,1 0,2 0,3 

fiuxon velocity, dimensionless units 



VIII. SUMMARY 

In summary, we analyzed the measurement process of 
a quantum state of a flux qubit by solitons propagating 
in an underdamped JTL. The coupling between the qubit 
and the JTL is inductive. We focused on the regime in 
which the information about the measured state is stored 
in the delay time of JTL solitons passing by (scattering 



FIG. 12: Different contributions to the error of a measurement 
for Ljtl = 25A.7, T = 20 mK, a = 0.001, 7c = 0.6 fiA and 
k = 0.05. The qubit stays at the symmetry point during the 
measurement. 



from) the qubit. If the qubit is in an energy eigenstate far 
from the symmetry point its persistent current induces 
an external magnetic flux in the JTL which serves as 
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a scattering potential for the JTL solitons. Two eigen- 
states induce different (opposite) scattering potentials, 
thus producing different delay times. A similar measure- 
ment mechanism works when the qubit is prepared at the 
symmetry point but, due to strong coupling, each passing 
soliton pushes the qubit far from the symmetry point. 

We have demonstrated that the delay times are longer 
for lower magnetic energies and velocities of the flux- 
ons and lesser dissipation in the JTL. Since the Joseph- 
son penetration depth is fixed by the design require- 
ment A « 2a the only possibility to reduce the mag- 
netic energy of the soliton is to decrease the Josephson 
junction critical currents Ic- For low critical currents 
Ic — 0.6 and Ic — 2.0 /iA the delay times are in the 
range td = 0.1 — lOw"^ which can be reliably detected by 
the RSFQ delay detector. The major sources of the mea- 
surement errors are the propagation time uncertainty of 
a fluxon due to thermal noise in the JTL and the intrin- 
sic qubit relaxation. Non-adiabatic transitions between 
the energy eigenstate induced by fluxons can cause an 
additional error if the qubit is measured at the symme- 
try point. For a JTL consisting of 50 elementary cells, 
dissipation strength a = 0.001, temperature T = 20 mK, 
fluxon velocity uq = 0.2, and coupling coefficient fc ~ 1 



the measurement errors of the qubit by a single fluxon 
are 9% and 26% for 1^ = 0.6 ^A and Ic = 2.0 /iA, re- 
spectively. In order to increase the signal-to-noise ratio 
one can use many fluxons which results in the improved 
accuracy of measurement exceeding 99% for the qubit 
far from the symmetry point and 90% at the symmetry 
point. 

For weak qubit-JTL coupling, k <^ 0.01, a qubit pre- 
pared at the symmetry point stays there all the time and 
induces no magnetic field in the JTL. In this case the 
measurement can be based on the fluxon scattered by 
the potential associated with the change of the effective 
inductance of that JTL cell which is coupled to the qubit 
loop. We found that a further reduction of the critical 
current of the Josephson junction is required for an ef- 
fective measurement in this regime. 



Acknowledgments 

We thank C. Hutter, A. Poenicke and M. Siegel for 
stimulating discussions and support. The work was sup- 
ported by the EU Specific Target Project RSFQubit. 



[1] D. V. Averin, K. Rabenstein, and V. K. Semenov, Phys. 

Rev. B 73, 094504 (2006). 
[2] V. K. Kaplunenko and A. V. Ustinov, Eur. Phys. J. B 

38, 3 (2004). 

[3] J. E. Mooij, T. P. Orlando, L. Levitov, L. Tian, C. H. 

van der Wal, and S. Lloyd, Science 285, 1036 (1999). 
[4] T. P. Orlando, J. E. Mooij, L. Tian, C. H. van der Wal, 

L. Levitov, S. Lloyd, and J. J. Mazo, Phys. Rev. B 60, 

15398 (1999). 

[5] K. K. Likharev, Dynamics of Josephson Junctions and 
Circuits. New York: Gordon and Breach (1996). 

[6] D. W. McLaughlin and A. C. Scott, Phys. Rev. A 18, 
1652 (1978). 

[7] A. V. Ustinov, Appl. Phys. Lett. 80, 3153 (2002). 
[8] B. A. Malomed and A. V. Ustinov, Phys. Rev. 69, 064502 
(2004). 

[9] Z. Hermon, E. Ben-Jacob, and G. Schon, Phys. Rev. B 

54, 1234-1245 (1996). 
[10] A.M. Savin, J. P. Pekola, T. Holmqvist, J. Hassel, L. 

Gronberg, P. Helisto, A. Kidiyarova-Shevchenko, Appl. 

Phys. Lett. 89, 133505 (2006). 

[11] J. Hassel, H. Seppa, and P. Helisto, cond-mat/05r0189 
[12] A. Lupascu, C.J. P.M. Harmans, and J.E. Mooij, Phys. 

Rev. B 71, 184506 (2005). 
[13] D.V. Averin in: "Quantum Noise in Mesoscopic 

P hysics" , E d. by Yu.V. Nazarov, (Kluwer, 2003) p. 229; 

cond-mat/0301524 
[14] Yu. Makhlin, G. Schon, and A. Shnirman, Physica C 368, 

276 (2002). 

[15] E. Joergensen, V.P. Koshelets, R. Monaco, J. Mygind, 



M. R. Samuelsen, and M. Salerno, Phys. Rev. Lett. 49, 
1093 (1982). 

[16] M. Salerno, E. Joergensen and M.R. Samuelsen, Phys. 

Rev. B 30, 2635 (1984). 
[17] N. Rosen and C. Zener, Phys. Rev. 40, 502 (1932). 
[18] I. Chiorescu, Y. Nakamura, C.J. P.M. Harmans, and 

J.E. Mooij, Science 299, 1869 (2003). 
[19] L. Gronberg, J. Hassel, P. Helisto, M. Kiviranta, 

H. Seppa, M. Kulawski, T. Riekkinen, and M. Ylilammi, 

Proceedings of the International Superconductive Elec- 
tronics Conference (2005). 
[20] P. Bunyk, A. Kidiyarova-Shevchenko, P. Litskevitch, 

IEEE Trans. Appl. Supercond. 7, 2697 (1997). 
[21] T. Ortlep and F. H. Uhlmann, to be published in IEEE. 

Tran. Applied. Supercond., (2006). 
[22] T. Flippov, M. Znosko, Proceedings of the International 

Superconductive Electronics Conference, (1999). 
[23] Q. Herr, D. Miller, J. Przybysz, Supercond. Science and 

Technology 19, S387 (2006). 
[24] S.B. Kaplan, A.F. Kirichenko, O.A. Mukhanov, S. Sar- 

wana, IEEE Trans. Appl. Supercond. 11, 978 (2001). 
[25] M. Terabe, A. Sekiya, T. Yamada and A. Fujimaki, to 

be published in IEEE Trans. Appl. Supercond. (2006). 
[26] P. M. Aziz, H. V. Sorensen, J. Van Der Spiegel, Signal 

Processing Magazine IEEE, 13, 61 (1996). 
[27] T. Ohki, A. Savin, J. Hassel, L. Gronberg and 

Anna Kidiyarova-Shevchenko, to be published in IEEE 

Trans. Appl. Supercond. (2006). 



